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Sound emergence in clarinetlike instruments is investigated in terms of instability of the static regime. 
Various models of reed-bore coupling are considered, from the pioneering work of Wilson and Beavers 
["Operating modes of the clarinet", J. Acoust. Soc. Am. 56, 653-658 (1974)] to more recent modeling 
including viscothermal bore losses and vena contract a at the reed inlet. The pressure threshold above 
which these models may oscillate as well as the frequency of oscillation at threshold are calculated. In 
addition to Wilson and Beavers' previous conclusions concerning the role of the reed damping in the 
selection of the register the instrument will play on, the influence of the reed motion induced flow is 
also emphasized, particularly its effect on playing frequencies, contributing to reduce discrepancies 
between Wilson and Beavers' experimental results and theory, despite discrepancies still remain 
concerning the pressure threshold. Finally, analytical approximations of the oscillating solution 
based on Fourier series expansion are obtained in the vicinity of the threshold of oscillation. This 
allows to emphasize the conditions which determine the nature of the bifurcation (direct or inverse) 
through which the note may emerge, with therefore important consequences on the musical playing 
performances. 



PACS numbers: 43.75.Pq [NHF] 



I. INTRODUCTION 

Sound production in reed wind musical instruments is 
the result of self-sustained reed oscillation. The mechan- 
ical oscillator, the reed, acts as a valve which modulates 
the air flow entering into the instrument, by opening and 
closing a narrow slit defined between the tip of the reed it- 
self and the lay of the mouthpiece. The phenomenon thus 
belongs to the class of flow-induced vibrations, which has 
been extensively studied both theoretically and experi- 
mentally (see, for example, Blevins^). Various regimes 
can occur in such systems: static regime, periodic oscil- 
lating regimes, and even complex chaotic behaviors. 

A first step to study this kind of oscillator, is to ana- 
lyze the stability of the trivial solution, the equilibrium 
position of the reed, in order to find a threshold of insta- 
bility associated with a set of control parameters defining 
the embouchure (reed, mouthpiece, and player) and the 
instrument itself. As an output of the threshold analysis, 
a mouth pressure threshold can be found. It is relevant 
for the musical playing performance: An estimation of 
the threshold value is a first evaluation of the ease of 
playing. In this paper the threshold of instability of the 
reed of a clarinetlike system is firstly investigated theo- 
retically by revisiting the important pioneering work of 
Wilson and Beavers^. Second, small oscillations beyond 
the instability threshold are considered. 

To answer the question whether the reed equilibrium 
is stable or not, several theoretical methods are available 
and have been used to study reed wind instruments. Us- 
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ing the feedback loop analogy, it is known as the free 
oscillation linear stability problem in a closed loop ob- 
tained when the nonlinear component of the loop is lin- 
earized around the trivial solution. Then the linear sta- 
bility of the reed can be studied with respect to each 
resonance of the input impedance of the resonator in fre- 
quency domain (see, for example, Refs.[l,[l,[5|,@,0 and H 
for reed instruments, Ref. |9| for lip-reed instruments and 
Ref . [l3 for a generic type of reeds) . Using the dynam- 
ical system representation where each resonance of the 
resonator is described as a simple second order oscillator 
in time domain, the linear stability analysis requires to 
solve eigenvalue problems and to analyze the sign of the 
eigenvalues real part (see, for example, Ref. 11 for lip- 
reed instruments and Ref. IH for vocal folds) . The the- 
oretical results can then be compared with experimental 
ones coming from artificial mouths by playingthem as 
gently as possible (see, for example, Ref s . Isl. Iisl and fl3 
for reed instruments; Ref . I 111 for lip reed instruments; 
and Refs. [13, [H, M and [ijfor vibrating vocal folds). 
The first attempt to derive theoretically the spectrum of 
reed instruments beyond the instability threshold is due 
to Worman^^. His results were at the origin of several 
works such as Refs. Ei, S [HI. 

Despite a rather simple description of physical phe- 
nomenon, pioneer theoretical and experimental results 
concerning mouth pressure and frequency at oscillation 
threshold were obtained by Wilson and Beavers^, show- 
ing the important role of reed damping determining the 
clarinetlike or lingual pipe organlike behaviour. 

The theory of Wilson and Beavers (hereafter denoted 
WB) is presented and discussed in Sec. [Til Then, model 
improvements from literature are added in the theory in 
order to try to reduce discrepancies between WB exper- 
iments and the theoretical results (Sec. IIIip . Section HVl 
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is devoted to a study of small oscillations beyond the 
instability threshold according to the direct bifurcation 
behavior of the clarinetlike instruments. Finally perspec- 
tives are discussed in Sec. [Vl 



II. WILSON AND BEAVERS THEORY 

The model and the principal results of WB are re- 
minded. Particular attention is brought to the numerical 
estimation of thresholds when the frequency approaches 
the reed resonance frequency, where some discrepancies 
with WB results are found. Finally the effect of the reed 
damping on playability is emphasized. 



A. Basic physical model 



bore representing a simplified clarinet body, an expres- 
sion given by Backus^: 

2Q 

(2) 

where p = —1. c, p, L, S and Zc = pc/S are the wave 
speed in free space, density of air, bore length, bore cross 
section, and characteristic impedance, respectively. The 
quality factor Q is assumed by WB to be frequency in- 
dependent, implying a damping proportional to the fre- 
quency. 

As it will be seen later, this assumption can be dis- 
cussed, and improved models will be used, because pres- 
sure thresholds are strongly infiuenced by bore losses 
(which are directly linked to the value of parameter Q). 



The physical model used by WB is reviewed with some 
comments related to more recent literature. It is based 
on the description of three separate elements: the reed, 
the bore, and the airfiow. The model used here is clas- 
sical and extremely simplified, but is also proven to be 
efficient in order to reproduce self sustained oscillations 
(for sound synthesis examples, see Refs. [22I and[23|). Ad- 
ditional elements will be discussed in the next section. 



1. Reed 

Based on the fact that reed displacement occurs in the 
vertical direction mainly without torsion, WB, among 
many authors, assumed a single degree of freedom mo- 
tion. Reed-lip-mouthpiece system is thus modeled as a 
lumped second-order mechanical oscillator with stiffness 
per unit area K, damping parameter g^, and natural an- 
gular frequency cOr, driven by the pressure drop Pm —p{t) 
across the reed, with an inward striking behavior: 

S + QrUJr^ + UJ^r (^(^) " ^o) = ^ {p{t) - Pm) , (1) 

p{t), Pm, y{t), and yo, being the mouthpiece pressure, 
the blowing pressure, the tip opening [denoted a{t) in 
WB's paper^], and the tip opening without any pressure 
difference, respectively. Pm is assumed to be constant. 

Avanzini et al^ numerically showed that this lumped 
model is valid for a small vibration theory where only the 
interaction between bore resonances and the first fiexu- 
ral mode of the reed is investigated. Measured transfer 
functions of a reed mounted on a mouthpiece also show 
a two degree of freedom response^^. 



3. Airflow 

As noted by Hirschberg^, in the case of clarinetlike 
instruments, the control of the volume fiow by the reed 
position is due to the existence of a turbulent jet. Indeed, 
a jet is supposed to form in the mouthpiece (pressure 
Pjet) after fiow separation from the walls at the end of 
the (very short) reed channel. Neglecting the velocity 
of air fiow in the mouth compared to the jet velocity 
'Ujet and assuming a downward air fiow {v^et > 0), the 
Bernoulli theorem applied between the mouth and the 
reed channel^^ leads to 

Pm = Piet ^ ^pvf^t' (3) 

Assuming a rectangular aperture of width W and height 
y{t), the volume fiow u across the reed channel can be 
expressed as follows: 

U{t) = Wy{t)^^^Pm-p;et(t). (4) 

Since the cross section of the mouthpiece is large com- 
pared to the cross section of the reed channel, it can be 
assumed that all the kinetic energy of the jet is dissi- 
pated through turbulence with no pressure recovery (like 
in the case of a free jet). Therefore, pressure in the 
jet is (assuming pressure continuity) the acoustic pres- 
sure p{t) imposed by the resonator response to the in- 
coming volume fiow u. This model is corroborated by 
experiments^^. Similar descriptions are used for double- 
reed instruments^^ and buzzing lips^. 



2. Bore 

The behavior of the acoustical resonator is determined 
by an input impedance relationship between acoustic 
quantities in the mouthpiece [acoustic pressure p(t) and 
volume fiow u{t)^ or P{(jo) and U{u)^ respectively, in 
the frequency domain]. WB assumed, for a cylindrical 



B. Characteristic equation and instability threshold 

The conditions for which self-sustained oscillations be- 
come possible are sought, that is, the minimum value of 
blowing pressure required for the static regime to be un- 
stable is investigated for a given configuration of the ex- 
periment {ur, Qrj and L being fixed). Common methods 
of linear stability analysis^^ are used in this study, and 
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solutions having time dependence exp {jcot) are sought. 
Cancehation of the imaginary part of uj corresponds to an 
oscihation that is neither damped nor amphfied: it char- 
acterizes the instabihty threshold of static regime. At- 
tention is drawn to the fact that this quantity may differ 
from the oscillation threshold depending on the nature 
of the bifurcation, studied in Sec. lIV.Bi As a language 
abuse oscillation threshold is often used instead of insta- 
bility threshold. 

Assuming small vibrations around equilibrium state 
(mean values of y and p are yo — Pm/K and 0, respec- 
tively), the volume flow relationship (|4]) is linearized. Di- 
mensionless quantities are introduced here: 6, 3^e, and D 
are the dimensionless frequency, input admittance and 
the reed transfer function, respectively. 



and D{e) 



1 



i + jqre-e-^' 

There are two dimensionless control parameters: 



7 = ^ and C = ZM\I^- 
Kyo V Kp 



(5) 



(6) 



7 is the ratio between mouth pressure and the pressure 
required to completely close the reed channel in static 
regime, while ( mainly depends on mouthpiece construc- 
tion and lip stress on the reed and is linked to the max- 
imum flow through the reed channel {( equals quantity 
2p in Ref. 2). 

Linearization of Eq. (|4]) leads to the so-called charac- 
teristic equation: 



which can be split into real and imaginary parts: 

Im(3^e(^)) =Cv^Im(D(^)) , 

1-7 



Re{ye{e)) =CV7 Re{D{0)) 



27 



(7) 

(8) 
(9) 



At last, a dimensionless length krL = cOrL/c is intro- 
duced. 



C. Numerical techniques 

The unknowns 0, e satisfying Eq. ([7j) are numer- 
ically determined for a range of bore lengths, parame- 
ters {qrX^^r) being set. They correspond to frequency 
and mouth pressure at instability threshold of the static 
regime. When various solutions exist for a given con- 
flguration due to the interaction of the reed resonance 
with the several bore resonances, the threshold observed 
experimentally by increasing the blowing pressure is the 
one having the minimum value of 7. 

The characteristic equation is transcendental and may 
have an inflnite number of solutions. Zero flnding is done 
using the Powell hybrid method^^, which combines the 
advantages of both Newton method and scaled gradient 
one. A continuation technique is adopted to provide an 



initial value to the algorithm: the flrst resolution is done 
for very high values of L {krL ^30), i.e., for a nearly 
non resonant reed (i.e., with neither mass nor damp- 
ing, hereafter denoted the massless reed case), where 
f {2n- l)c/4L and 7 1/3 (with n G N*). Bore 
length is then progressively decreased and zero flnding 
for a given value of L is initialized with the pair (^, 7) 
solution of the previous solving (bore slightly longer). 
Depending on the initialization of the flrst resolution 
{krL = 30), it is possible to explore the branches asso- 
ciated with the successive resonances of the bore. When 
reed resonance and bore antiresonance get closer to each 
other {krL n7r{n G N)), fast variation of the pressure 
threshold requires adjustment of the step size. 



D. Results 

Two kinds of behavior can be distinguished. For 
strongly damped reeds (Fig.[T]), the threshold frequencies 
(dashed lines, result from Ref. ^ always lie either near 
the reed resonance (^ = 1) or near the flrst impedance 
peak frequency of the pipe (hyperbola = tt/ {2krL), not 
represented in the flgure for readability), corresponding 
to the flrst register of the instrument. When the length 
L decreases, pressure threshold gradually reduces from 
values assumed for the massless reed model to a mini- 
mum point for krL ~ 7r/2, and then strongly increases 
as the pipe becomes shorter. When increasing 7 from 
0, the loss of stability of the static regime may give rise 
to an oscillating solution which always corresponds to 
the flrst register since the instability thresholds of the 
higher registers occur for higher values of mouth pres- 
sure. On the contrary, considering now slightly damped 
reeds (Fig. [2]), emerging oscillations can occur near higher 
pipe resonances. Indeed, for certain ranges of L, the 
pressure threshold associated with one particular higher- 
order register is lower than the pressure required to drive 
the air column in the flrst register. This lowest pres- 
sure threshold is associated with the acoustic mode, the 
natural frequency of which being the closest to the reed 
resonance. 

These results show the influence of reed damping on 
selecting 

• a clarinetlike behavior with heavily damped cane 
reed: preference is given to the chalumeau register, 
i.e., the lowest register, or 

• a lingual organ pipe behavior with very lightly 
damped metallic reed: tuning is performed by ad- 
justing reed vibrating length by means of a wire 
adjusted on reed, see Ref. [32| for further details. 

This noticeable conclusion of Wilson and Beavers" paper 
gainsays von Helmholtz^^ , who stated that the differences 
of behavior are linked to the mass of the reed. 



E. Discussion 

Despite the important results obtained by WB, dis- 
crepancies concerning numerical results and model limits 
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FIG. 1. Dimensionless threshold frequencies (top) and pres- 
sure (bottom) for a strongly damped reed: qr = 0.4, fr = 
750 Hz, p = 0.065. Results (dashed lines) and measurements 
(squares) from Ref . Q; our numerical results (solid lines). 



earization beyond model limits. 

The fact that the dimensionless mouth pressure 7 re- 
mains lower than 1 implies the existence of a maximum 
threshold frequency for the oscillation. We assume that 
the maximum value ^max 

is associated with the maximum 
value of 7 (i.e.,1). No theoretical proof is given here, but 
this assumption is observed in all numerical results ob- 
tained. Then the solution ^max is given by 

:^^e(^max) = C-D(^max), (10) 

which leads to relation 

= CRe(2e(^max)) >0 (11) 



1 



92 

^max 



as the bore is a passive system. As a conclusion, it is not 
possible to play sharper than a frequency slightly flat- 
ter than the reed-lip-mouthpiece system resonance fre- 
quency, the deviation to this latter being linked to the 
losses in the bore. 
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FIG. 2. Dimensionless threshold frequencies and pressure 
for a lightly damped reed: qr = 0.008, fr = 700 Hz, p> = 
0.05. Results (dashed lines) and measurements (squares) from 
Ref. Q; our numerical results (solid lines). 



can be pointed out. In a validation phase of our numer- 
ical algorithms, our results (see Fig. [T] and [2]) were com- 
pared to their ones (the data being extracted from their 
article), solving exactly the same equation. Differences 
appear between numerical results concerning threshold 
frequencies, approaching reed resonance when the length 
decreases. Investigations pointed out model limits the 
authors do not take into account: linearization of flow 
relationship is valid only while reed channel is not closed 
at rest, i.e., when < Ky^. Using the linear form 
for higher values of Pm would be meaningless, even for 
free reed aerophones: the opening function (linked to 
the reed displacement) taking part in flow calculation 
can never be negative. For instruments for which reed 
beats against the mouthpiece, reed channel is completely 
closed and then sustained oscillations cannot occur for 
bore lengths where WB theory predicts pressure thresh- 
old above static beating reed pressure by extending lin- 



F. Minimum pressure threshold: Improved playability for 
interacting resonances 

For each value of ^r, there exists one or more ranges of 
bore lengths where playability is greatly improved. In- 
deed, pressure threshold curves show a minimum for a 
certain value of /c^L, denoting an increased easiness to 
produce the note corresponding to this length. Associ- 
ating a clarinet mouthpiece with a trombone slide, in- 
formal experiments conflrm that it is easier to produce 
some notes than other ones. Analytical approximated ex- 
pressions for the minimum have been investigated. Un- 
der the assumption that this minimal value is obtained 
for an emerging frequency located close to the reed reso- 
nance, and therefore is mainly determined by reed damp- 
ing, bore losses can be ignored Re(3^e(^)) = 0, Eq. (|9]) 
leading thus to: 



7 



1 + 2Re (L>(i9)) * 



(12) 



In the massless reed model [p{p) = 1), threshold pres- 
sure is equal to 1/3 and frequencies correspond to fre- 
quencies for which the imaginary part of the bore input 
impedance vanishes, which is consistent with results al- 
ready published^^. The minimum pressure threshold oc- 
curs at a maximum of Re {D{0)) : 



Re{D{0)) 



1 



(13) 



obtained for = ^1 — qr (which is consistent with the 
approximation 6> :^ 1), thus, 

qr{2 - qr) 



70 



2^qr{2-qr)' 



(14) 



nearly proportional to qr for lightly damped reeds. 

In order to evaluate the effect of bore losses on this 
minimum pressure threshold, a one-mode resonator with 
losses is now considered: 



ye{oo)=Yn(l^jQn 



CO 



CO 



(15) 
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where Yn is the minimum amphtude of the admittance 
and Qn is the quahty factor; Eqs. (|8]) and (|9]) become 



7 



1 



2V7 '^'(1- 



- ^2)2 + ^q^gy 



(16) 
(17) 



Ignoring the variation of Yn with the length of the bore, 
the derivation of Eq. ([16]) with respect to krL leads to 
a minimum value of the function 7 = fik^L) for 9'^^^ = 
I — qr and 7niin solution of equation 



Yn 



1 



■7n 



1 



CV7min 27min <7r(2 - ^r) ' 

which the first-order solution is given by 



70 



1 



obtained for 



1 _ *: 



C 



2Q. 



(18) 



(19) 



(20) 



For an open/closed cylinder uj^ 
result is 



ri 



(2n-l)- 1 + 



2 



= (2n - l)7rc/2L, the 

1 c 



2Qn '^YnQn's/Or 



(21) 

Typical values Fn = 1/25, C = 0-4, and = 0.4 lead to 
an increase in 7min relative to 70 of about 8%, confirming 
the preponderant effect of reed damping on the minimum 
pressure threshold. 

In order to understand how coupling acoustical and 
mechanical resonances could reduce pressure threshold 
over a wider bore length range, the neighborhood of the 
previously mentioned minimum has been studied. Using 
again a single acoustical mode for the calculation, deriva- 
tion of a parabolic approximation was possible. Writing 

7 = 7min(l+S:^), = Cin+^ ^ud UJn = (c^n)min (1 + , 

with £, 6, and v small quantities, the Taylor expansion 
of Eqs. ([T6|) and ([T7|) with respect to these values leads 
to the next relationships: 



2qlY„Qn{5 - 2v) = -5Cy 



(22) 
(23) 



Finally, near the minimum pressure threshold, depen- 
dence on bore length is given by 



7 = 7n 



2<7r 



,3/2 . 



krLiYiin 



(24) 

In a first approximation, the aperture of the approxi- 
mated parabola, thus the width of the range for which 
oscillation threshold is lowered, is mainly controlled by 



the musician embouchure, i.e., by reed damping and lip 
stress on the reed. This means that, thanks to its em- 
bouchure, the player can expect an easier production of 
tones for certain notes. 

For a lossy cylindrical open/closed bore, modal ex- 
pansion of input impedance gives YnQn = uJnL/2c = 
(2n — l)7r/4 so that bore losses do not seem to have a 
great infiuence on playing facility, at least when consid- 
ering minimal blowing pressure 7min {Qn does not appear 
alone in first order calculation). On the contrary, they 
are essential for understanding the extinction threshold 
phenomenon^^, i.e., when the reed is held motionlessly 
against the lay. 



III. MODEL IMPROVEMENTS 

Last four decades have been fruitful in physical mod- 
eling of musical instruments, especially for single reed 
instruments. Pipes have been the focus of a great num- 
ber of studies since Benade^, as well as the description of 
peculiarities of the fiow (Backus^, Hirschberg^, and Dal- 
mont et al!^^). The aim here is to try to reduce discrepan- 
cies between WB experiments and theory, based on some 
of those investigations which look relevant to the study of 
oscillation threshold. In a real clarinet-player system, the 
vocal tract may have an effect on the oscillation threshold 
(see, e.g., Ref . IstI) . However, the use of porous material 
in the WB apparatus suggests that acoustic interaction 
between the chamber and the mouthpiece has been can- 
celed or at least strongly reduced. Therefore we do not 
discuss the effect of the upstream part. 



A. Viscothermal losses model and vena contracta 

It should be noticed that there is a gap in pressure 
threshold values between experiment and theory in WB 
article. This occurs even for long bores when reed dy- 
namical behavior should not deviate from the ideal spring 
model (because of an emerging frequency much smaller 
than uor). For that case, Kergomard et al^ provided 
an approximated formula taking into account both reed 
dynamics and bore losses that can be extended to 



7 ■ 



el , 2Re(3^e(^n)) 



3-€ 



3V3C 



(2n-l) 



'2krL 



■ (25) 



On corresponding to the nth resonance frequency of the 
bore. In comparison with ideal model (lossless bore and 
massless reed, i.e., 7 = 1/3), additional corrective terms 
are considered in Eq. ([25|) , one lowering pressure thresh- 
old due to the collaboration of the resonant reed, the 
other one requiring higher blowing pressure due to dissi- 
pation in the bore. According to this approximated ex- 
pression, pressure thresholds depend on the mouthpiece 
parameter and on bore dissipation at playing frequency. 

Now focus is done on using realistic values of ( and 
ye{0), assuming acoustic losses in clarinetlike bore to be 
due mainly to viscothermal dissipation (Ref. 36). Other 
kinds of losses such as nonlinearity localized at the open 
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end of a tube are negligible since study is done at oscilla- 
tion threshold, i.e., for very small amplitude oscillations. 
Simpler model ("Raman's model") has been recently in- 
vestigated by Dalmont et al^ and led to satisfactory 
oscillation threshold of the fundamental register of the 
clarinet when operating frequencies are much lower than 
the reed natural frequency (the reed being considered as 
an ideal spring). In the present study, the magnitudes of 
the higher order impedance peaks need to be correctly es- 
timated so that neither the loss model of Raman nor that 
of Backus^ may be realistic enough. Pressure thresholds 
would be inaccurate when higher-order modes of bore 
oscillate first. 

The standard formula for the input impedance will be 
hence considered: 



Ze(u) = jZctdiiilkL) with jk = — ^r]\ — (26) 

c \ c 



where 77 is a coefficient quantifying the viscothermal 
boundary layers, equals to 0.0421 in mks units for a 
7 mm radius cylinder. This model introduces dissipa- 
tion and dispersion and leads to a zero valued impedance 
at zero frequency, which is still consistent with the lin- 
earization of ffow relationship with a zero mean value of 
acoustic pressure. In this expression, viscothermal effects 
are ignored in the characteristic impedance (see Ref . 34) . 
Equation ([26|) leads to peak magnitudes decreasing with 
length L, whereas they are not sensitive to L in Backus 
empirical expression. A direct consequence is that the 
pressure threshold increases as the resonator lengthens. 

Another effect may occur and modify pressure thresh- 
old. Hirschberg^ and Hirschberg et al^^ brought to at- 
tention on vena contracta phenomenon: due to sharpness 
of edges, ffow separation may result in the formation of a 
free jet in reed channel, this contraction effect resulting 
in a jet cross-sectional area smaller than the reed chan- 
nel opening. Recent investigations^^ applying the lattice 
Boltzmann method to the reed channel confirm previous 
experiments^. The assumption of constant vena con- 
tracta is valid in some specific cases, for short channel 
geometry, but there is no proof yet that this effect is sig- 
nificant for a real clarinet mouthpiece. Here small vibra- 
tions of the reed near oscillation threshold are considered, 
which may induce little infiuence of the fiow unsteadiness 
on the measured pressure threshold. So it is possible^^ to 
include vena contracta phenomenon by reducing the area 
of the reed channel Wy{t) by a coefficient approximately 
0.6, i.e., by multiplying ( by this coefficient. 

Numerical investigations point out (see Fig. [3]) that 
taking into account realistic losses and vena contracta 
phenomenon reduces discrepancies in pressure thresh- 
old, especially for high values of bore length and for 
a strongly damped reed, i.e., when reed dynamics has 
a small inffuence on oscillation threshold. Nevertheless 
pressure values experimentally obtained by WB still re- 
main quite higher than the ones corresponding to the 
modified model, while frequencies are unaltered. 

An attempt to explain the discrepancies in threshold 
measurements will be provided in Sec. HVl 




FIG. 3. Results with the model considering vena contracta 
phenomenon (thin line) and viscothermal losses (thick line) 
(same conditions as in Fig. [T]). WB results are also recalled 
(dashed line). 



B. Reed motion induced flow 

The infiuence of the reed is not limited to its resonance. 
Since Nederveen^^ and Thompson^^, it is proved that the 
fiow entering through the reed opening is divided in one 
part exciting the resonator and another part induced by 
the reed motion. In fact, the vibration of the surface of 
the reed produces an additional oscillating fiow. Thus 
the entering fiow U can be written as 



(27) 



where Sr is the effective area of the vibrating reed related 
to the tip displacement y{t). Alternately, a length Al 
can be associated with the fictitious volume where the 
reed swings. Dalmont et al}^ reported typical values 
of 10 mm for a clarinet. Nederveen^^ linked A/ to reed 
strength (or hardness): Al may approximately vary from 
6mm (strong reeds) to 9 mm (softer reeds). These values 
being small compared to clarinet dimensions, the reed 
motion induced fiow can be considered through a mere 
length correction in common work, but its infiuence on 
the interaction between acoustic resonator and reed is 
not negligible on the threshold frequency, as it is studied 
now. 

In a first step, this effect is considered separately, all 
losses being ignored {qr = and 7^ = 0). Equation ([27|) 
coupled to Eq. ([7j) leads to the following system: 



- Im(3^e) = krAl 



0^ 



1-7 ^ 1 - 6>2 



(28) 
(29) 



As seen previously, several frequency solutions exist for 
different pressure thresholds 7. Examination of Eq. ([29|) 
as a function 7 = f{6) (for ^ < 1) reveals that the so- 
lution having the lowest threshold is the one being the 
closest to the reed frequency. 

Approximations can be derived in some situations. 
When playing close to a bore resonance frequency On <C 
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0.0 



FIG. 4. Comparing reed motion induced flow effect and 
reed damping effect (A/ = 12 mm and other conditions as 
in Fig. [T]): reed motion induced flow only (dashed lines), 
reed damping only (dash-dotted lines), and both effects (plain 
lines). WB experimental results are recalled (squares). 



1, the right-hand term in Eq. (|28|) is small, so that, with 
3^e = —j cot {6krL), the reed motion induced flow acts 
merely as a length correction: 



Al 



1 



where On = 



(2n- l)7r 



(30) 



This approximation is valid when the considered bore 
frequency On remains smaller than unity. An equivalent 
approximated length correction can be derived for the 
effect of reed damping on frequencies from Eq. (|8]): 



Ala 



C 



(31) 



A trial and error procedure has been performed to adjust 
the oscillation frequency. Very few iterations were needed 
to exhibit a value for the reed motion induced flow (A/ 
12mm) higher than for the reed damping {Alq c:^ 2mm) 
in the conditions of Figs. [T] and [H When acoustical and 
mechanical resonances are very close {0 = 1 — e and On = 
1 — en)j a second-order expression can be deduced: 



(32) 



the apparition of a square root being typical of mode cou- 
pling, making the achievement of analytical expressions 
difficult. Then, when the bore length decreases enough 
so that one of its resonances increases above the reed one 
{On > 1), the oscillation frequency approaches the reed 
one until the intersection point disappears for krL = utt 
(see Fig. [5]). Near the reed resonance, first-order approx- 
imations can be derived: 




c^l — -krAl tan (krL)^ 
7 ^krAl tan {krL). 



(33) 
(34) 



FIG. 5. Graphical representation of Eq. (|28|) giving oscilla- 
tion frequency at threshold: left-hand term — Im {y) (thick 
lines), right-hand term krAW/{l — 6^) (thin line), and solu- 
tions (markers). 



These expressions are valid if tan (/c^L) > 0, i.e., krL > 
nn. According to Eq. ([29|) . when oscillation frequency ap- 
proaches fr, the pressure threshold decreases to zero con- 
trary to what happens when considering the reed damp- 
ing effect. 

Figures [4] and [6] show a numerical comparison of the 
respective effects of reed motion induced flow and reed 
damping. The reed motion induced flow with Al = 
12 mm in Fig. [4] and 5 mm in Fig. [6] adjusts the frequency 
deviation for both heavily and slightly damped reeds, 
even when approaching reed resonance, and is prepon- 
derant compared to the damping effect. An example of 
a mistuned length correction parameter {Al = 12 mm) 
is exhibited in Fig. [6l revealing a mismatch of oscilla- 
tion frequency. This latter seems to be mainly controlled 
by the reed motion induced flow effect when approach- 
ing the reed resonance frequency. Classical approaches 
(see, e.g., Ref. [4lh for the calculation of playing frequen- 
cies, ignoring the flow due to pressure drop, and search- 
ing for eigenfrequencies of the passive system including 
the bore and the reed only are thus justified. It can be 
noticed that Eq. ([28|) was already given by Weber^^ in 
the early 19th century (see page 216), assuming a reed 
area equal to the cylindrical tube section. This theory, 
used by several authors (see, e.g., Ref. 32), was discussed 
by von Helmholtz^^ and Bouasse^^, especially concerning 
the lack of explanation on the production of self-sustained 
oscillations. On the contrary, threshold pressure curves 
exhibit that both damping and additional flow have in- 
fluence on the pressure required for the reed to oscillate. 
So a combination of the phenomena has to be taken into 
account, none of them being negligible in the considered 
domain. 



IV. GOING BEYOND INSTABILITY THRESHOLD 

A. Linear stability analysis with modal decomposition 

Analysis of the instability threshold can be performed 
using complex frequency formalism. For a given config- 
uration of the whole system hore-reed-musician (L, 
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FIG. 7. Evolution of the complex eigenfrequencies as a func- 
tion of mouth pressure 7. r = 7 mm, cjr = Stt x lOOOrad/s, 
0.3, L = 16cm {Oi 0.53 and KL 2.96), and C 0.2. 



FIG. 6. Comparing reed motion induced flow effect and reed 
damping effect (conditions as in Fig. [2]): reed motion induced 
flow only (A/ = 5 mm: plain thin lines. A/ = 12 mm: dotted 
line in upper graph only), reed damping only (dashed thin 
lines), and both effects (plain thick lines). WB experimental 
results are reacalled (squares). 
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cjr, <7r, 7, and ^ being set), its complex eigenfrequencies 
Sn = j^n — <^n Can be determined. The imaginary part 
of Sn corresponds to the frequency, the real part an be- 
ing the damping of this mode, for the coupled linearized 
system close to the static equilibrium state. Classically, 
when the mouth pressure is below the oscillation thresh- 
old, all eigendampings an are positive, the static regime 
being stable. An oscillation may appear when at least 
one mode of the whole system becomes unstable, i.e., 
when at least one of the an becomes negative. Looking 
for instability threshold can be done by varying a bifurca- 
tion parameter (either L, 7, or () and examining the real 
part of computed eigenfrequencies. Two examples are 
shown in Figs. and [3 It is noticeable that the frequen- 
cies evolve only slightly with the bifurcation parameter 7 
and are close to eigenfrequencies of either the bore or the 
reed ( Im (jcj/cj^) 1). For the first example, the static 
regime becomes unstable for 7 ^ 0.28 and a frequency 
near the first resonance of the bore (dot-dashed curves). 
Other acoustic resonances (their frequencies being odd 
multiples of the first one) and reed resonance have higher 
oscillation threshold and remain damped for this configu- 
ration. For a longer tube (Fig. [8]), instability appears for 
7 :^ 0.3 at a frequency located near the third resonance 
of the bore [dotted line at Im {juo/uor) ::^0.8:^3x0.27], 
the first resonance becoming unstable for a larger mouth 
pressure. The frequency close to the reed resonance (solid 
curves) still remains damped. 

Calculations may be simplified and accelerated by us- 
ing a modal decomposition of the bore impedance Ze{uj): 
this allows for the characteristic equation to be written 
as a polynomial expression of jco, and optimized algo- 
rithms for polynomial root finding can be used. Modal 
expansion considers the N first acoustic resonances of the 



FIG. 8. Evolution of the complex eigenfrequencies as a func- 
tion of 7. r = 7mm, ujr = 27r x lOOOrad/s, qr = 0.3, 
L = 32cm {Oi 0.26 and krL 5.91), and C 0.2. 



bore: 

Zc 



- j tan 



c 



■ja{u)L 



(35) 
(36) 



where modal coefficients ujn and Qn can be deduced ei- 
ther from measured input impedance or from analytical 
expression [Eq. (|35|) ]. assuming a{u;) to be a slowly vary- 
ing function of the frequency. 

Comparison between direct calculation of oscillation 
threshold 7th using WB method and estimation using 
modal decomposition and complex eigenfrequencies com- 
puting has been done. Whereas the difficulties for the 
first method arise due to the transcendental characteris- 
tic equation, the second one requires calculation of eigen- 
values for various values of the mouth pressure 7, using 
an iterative search of the instability threshold. The num- 
ber of modes taken into account has been chosen such 
that resolution of 7th is less than 0.01, which is also the 
tolerance used for the iterative search. An example is 
given in Table [J for a heavily damped and strong reed. 
It shows very good agreement between the two methods 
for both 7 and 6. This validation allows the use of the 
complex frequency approach, which results in an efficient 
algorithm that can be easily applied to more complex res- 
onators whenever a modal description is available. 
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OcF OwB 


AO 






A7 


8.5 


u.xou u.xou 




0.43 


0.43 


1% 


2 


n 747 n 747 


0.1% 


0.30 


0.30 


0.6% 


1 


1.022 1.024 


0.2% 


3.86 


3.82 


1.0% 


0.81 


1.033 1.034 


0.1% 


8.77 


8.71 


0.7% 



TABLE I. Comparison of pressure threshold and oscillation 
frequency calculated using complex frequency formalism (in- 
dexed by CF) and Wilson & Beavers method (WB) for 
r — 7mm, ujr ^ 27i x 750rad/s, qr — 0.4 and C — 0.13. 



The technique is the harmonic balance applied to oscil- 
lations of small amplitudes. Calculations are done here- 
after by taking into account the reed dynamics in the 
volume flow relationship, which does not appear in the 
mentioned paper. Fourier series of the volume flow de- 
pends on Fourier components of signal Pm — p{t) and 
y{t). Assuming steady state oscillations with angular fre- 
quency CO, the dimensionless signals are written as 



Pit) 



U{t) = Uo^^ynPn 



(38) 



The writing of the characteristic equation for a sin- 
gle acoustic mode exhibits the behavior of the coupled 
oscillators: 

2 ■ f c , 1 - 7 

r 2 . 21 . 2c 2 ^ . AA 

' (37) 

Coupling realized by the flow in the reed channel mod- 
ifles the damping of the acoustic mode: in addition to 
the usual term (corresponding to viscothermal losses and 
eventually radiation) , damping is increased by a quantity 
related to mouth pressure and stress on the reed. This 
may be regarded to as a resistive acoustic behavior at the 
bore entrance. 

Assuming that a linearized model is still relevant dur- 
ing the growth of oscillations (before the amplitude sat- 
uration mechanism appears), this approach can be ex- 
tended to investigate the transient response of the cou- 
pled system. Characterizing the degree of instability 
of the system by a = min^ o^n, the slope of the curve 
(7 = /(7) gives an information about the instability de- 
gree of the system when the mouth pressure is slightly 
higher than the oscillation threshold. A great slope 
would correspond to a very unstable conflguration and 
a quick growth of oscillation, whereas nearly constant 
curve would lead to a small ampliflcation coeflicient and 
slowly rising vibrations and then to longer transient at- 
tack before stabilization of the magnitude of oscillations. 
Links between the computed eigenfrequencies of the cou- 
pled system presented here and the transient behavior 
have still to be investigated. 



njojt 



(39) 



where J'n = ye{nix}) and Dn = D{nu;) are the values of 
dimensionless bore admittance and reed transfer function 
for angular frequency nuj. The volume flow relationship 
is rewritten as 



uHt) = e {y{t)/yoY{i-p{t)). 



(40) 



Sustained oscillations of very small amplitude are stud- 
ied, assuming that pi is a non vanishing coeflicient con- 
sidered as a flrst-order quantity. Notations CE^ and 
Y(n\ are introduced: 



CE„ = yn/{c,^) + ^ - -D„, 



27 



(41) 



"^(^) ^ "^(m) ^ 



7 



(42) 



Cancellation of CE^ for given uo and 7 means that the 
characteristic equation ([7j) is solved for nuo and 7. Ex- 
panding Eq. (|4Q|) leads to 







2(1-7)E 



u^yn 



C27(l - 7) 

,{n-\-m)jujt 



Dn 



1-7 

27 



njujt 



,Pm^ 



iPnPmPq^ 



{n-\-m-\-q)jujt 



(43) 



B. Mouth pressure required to obtain a given (small) 
amplitude 

The previous sections of the paper deal with the stabil- 
ity of the static regime, looking for the condition to make 
a bifurcation possible. Some developments concerning 
the existence of oscillating regime above the threshold 
are derived now. Neither stability of oscillation nor tone 
deviation issue will be discussed here. 

Grand et al!^^ suggested the introduction of the lim- 
ited Fourier series of pressure in the massless reed case. 



It is here assumed that Pn is of order \n\ (with p-n = Pn) 
(see, e.g., Ref. 20). The continuous component of the 
volume flow is calculated up to order 2: 



^=(l-7)^ + EF(+n)K 



n/0 

PnPmPn+m 

(l-7)' + 2F.+n|pi|2 + o(p2) 



(44) 
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so that, considering uq to be real 

\P 



(1-7)2 



+ o{pi). (45) 



From Eq. (pS)) . the frequency (Nuj) is extracted for 
AT > 1: 



= 2(1-7) 



D 



N 



1-7 



Uoy-n 



27 ("7(1- 7) J 



PN 



1 

7 



iPnPmPN- 



(46) 



For N > 2, Taylor series expansion up to order N is 
applied: in the first sum, only terms corresponding to 
< n < N contribute at order N, while, in the second 
one, the terms to consider are the ones for which < 
n < N and < m < N — n. The component pn can be 
deduced from the sequence {pn)o<n<N- 



Pn 



1 



2(1-7) CE^ 



+ o(pf). (47) 



_ 0<n<Ar 

mPnPmPN —n—m 



0<n<N 
0<m<N-n 



As an example, for N = 2, the second sum being empty: 



P2 



2(1-7) CE2 



oipi). 



(48) 



As expected and in agreement with the so-called "Wor- 
man rule"^^, higher components appear to be higher or- 
der quantities: order 2 for p2, 3 for ps, 4 for p4, etc. 

Focus is now given to fundamental frequency. Calcu- 
lations are done up to order 3: 



0=2(1-7) 



Di 



1-7 



uoyi 



27 ejii - 7) 

+ E ^(i!!„)P«^'i-« 

E DnDmPnPmPl-n- 

1-7 



pi 



7 

.2(1-7) 



(49) 



upyi 

27 ("7(1 - 7) 
1 



pi 



+ 2F(_2^)P2Pt - -DiPi\pi\'{D, + 2DI) 



Replacing uq and p2 gives 



\P1? 



2(l-7)'CEi 



1-7 



(D? + 2|Z)i|2) 
(50) 



The limitation to the first harmonic method ignoring the 
infiuence of Pn>2 on the amplitude \pi \ leads to 



\Pi? 



2(l-7)2CEi 



23^1 p 



1 



-(D2 + 2|Di|2) 



(51) 



In their paper. Grand et al!^^ stated that the procedure 
used in a first simplified case can be applied to a model 
including reed dynamics provided that Zq{(jo) is replaced 
by Ze{u)D{u). Conclusion is not so straightforward as 
reed dynamics also interferes with the volume fiow rela- 
tionship contributing to a more complex expression for 
the first component amplitude. According to their work, 
Dn terms would only occur with in the expression of 
which is not the case. 
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FIG. 9. Comparison of computed with the harmonic bal- 
ance for small oscillations for various values of krL\ from 
thin to thick lines, KL = 9, 5, 4, 3, 2, 1.8, 1.5, 1.25 (A/ = 0, 
qr = 0.4, and fr = 1050 Hz). 



In Fig.[9]are shown the bifurcation diagrams for a heav- 
ily damped reed. From the fiattest tone {krL = 9.00) to 
the sharpest computed {krL = 1.25), oscillations with 
very small amplitudes are possible as the diagrams show 
Hopf bifurcations that are supercritical, i.e., direct, for all 
the computed cases. On the other hand, for cases where 
oscillation frequency is closed to the reed one (Fig. [TQ|) . 
a subcritical Hopf bifurcation occurs for a range of bore 
length, with a very small pressure threshold linked to 
the interaction of reed resonance with the second bore 
resonance: the bifurcation is there inverse. 

In accordance with Grand et a/.^^, it appears that the 
bifurcation is not always direct. There exist configura- 
tions for which computed bifurcation diagram shows a 
subcritical pitchfork, i.e., small oscillations for values of 
pressure below static regime instability threshold. The 
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FIG. 10. Comparison of \pi\ computed with the harmonic 
balance for small oscillations for 14 values of krL regularly 
spaced between 3.6 and 4.20 with a lightly damped reed (A/ = 
0, qr = 0.01, and fr = 1050 Hz). 



boundary between the two cases is not trivial to ex- 
plore analytically, requiring the mathematical study of 
Eq. ([5Q|) . Alternately numerical exploration of param- 
eter space can lead to some partial observations. One 
of them is that the bifurcation seems to be direct when 
the oscillation frequency is close to one of the bore reso- 
nances. 

The limit of calculations derived here needs to be em- 
phasized. When \pi\ tends to zero, the oscillation fre- 
quency is the one for which the characteristic equation 
is solved ( CEi = 0). Higher components (the second 
one here at the first-order approximation) have a non- 
negligible infiuence only if CE2 becomes small too when 
approaching oscillation threshold, i.e., if the character- 
istic equation accepts the solutions co and 2a; for two 
close values of mouth pressure: CEi and CE2 are small 
simultaneously. Inverse bifurcation may occur in the de- 
generate case of simultaneously destabilization of static 
regime for a frequency and its octave. This considera- 
tion can be extended to the more general case of CEi 
and CEat {N > 2) canceling an equal mouth pressure 
value. 

As a conclusion, the nature of the bifurcation depends 
on the roots of the characteristic equation ([7j): this gen- 
eralizes the results of Grand et al^. 



V. CONCLUSION AND PERSPECTIVES 

Two components of the volume velocity at the input 
of the resonator act on the oscillation threshold values: 



the first one is due to the pressure drop between mouth 
and mouthpiece, while the second one is due to the reed 
movement. Roughly speaking, the first one has a mainly 
resistive effect, either passive or active, the second one 
having a mainly reactive effect. This remark can be re- 
lated to the behavior of the threshold pressures and fre- 
quencies. Concerning the frequencies, the second effect 
is preponderant, at least for the cases studied by Wilson 
and Beavers, and it can justify the historical method due 
to Weber^^ regarding the playing frequencies as eigenfre- 
quencies of a passive resonator. Concerning the pressure 
thresholds, the fiow due to the pressure drop is essen- 
tial, and the model based on Bernoulli equation used by 
Wilson and Beavers^ seems to be satisfactory, but taking 
into account the fiow due to reed movement is necessary, 
and improves the results of the authors, mainly for lightly 
damped reeds. 

Discrepancies between Wilson and Beavers' experi- 
mental results and numerical ones still remain. Can this 
be due to the nature of the bifurcation at threshold? For 
a clarinetlike functioning, i.e., for strongly damped reeds, 
numerical calculation of small oscillations proves that it is 
supercritical, confirming the works ignoring reed dynam- 
ics, while it is not sure for lightly damped reeds. Further 
experimental investigations are planned to conclude on 
the nature of the bifurcation on real clarinets. Another 
unknown concerns the acoustic impedance of the air sup- 
ply device and of the mouth cavity. 

Concerning the possible generalization of this work, 
other kinds of either resonators or reeds should be stud- 
ied, as, e.g., in Ref. |8|. Some authors of the present paper 
recently obtained simplified theoretical results for an out- 
ward striking reed by using a similar method of investi- 
gation based on a single degree of freedom reed model^. 
Calculations with several DOFs for the reed remain to 
do, especially for lip reed instruments (see Ref. fllh . even 
if recent works dedicated to vocal folds have been done 
(see Ref. |l2|) . Finally, the study of transient sounds can 
be made easier with a good knowledge of the linearized 
functioning of reed instruments. 
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